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We develop a theoretical model to describe the dissipative dynamics of singlet-triplet (S-Tq) qubits in GaAs 
quantum dots. Using the concurrence experimentally obtained as a guide, we show that each logical qubit is 
coupled to its own environment because the decoherence effect can be described by independent dephasing 
channels. Given the correct description of the environment, we study the dynamics of concurrence as a function 
of the temperature, the constant coupling between the system and the environment, the preparation time, and 
the exchange coupling. Although the reduction of the environment coupling constant modifies the entanglement 
dynamics, we demonstrate that temperature emerges as a crucial variable and a variation of millikelvins signif¬ 
icantly modifies the generation of entangled states. Furthermore, we show that the exchange coupling together 
with the preparation time strongly affects the entanglement dissipative dynamics. 
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I. INTRODUCTION 

The development of quantum information processing have 
discovered different new techniques which are paving the way 
to accomplish quantum technological advancement.— Among 
these advances, spin qubits in quantum dots (QDs) is cer¬ 
tainly one of the most striking systems because of their po¬ 
tential scalability and miniaturization—^ Furthermore, elec¬ 
trical readout and control of single spins in quantum dots 
(QDs) have been proven very challenging, where spin block¬ 
ade and charge sensors enable the observation of single/two- 
spin dynamics 

More recently, a significant progress in implementing an al¬ 
ternative scheme in double quantum dots (DQDs) has been at¬ 
tained.— In this approach, singlet and triplet states of two elec¬ 
trons are used to represent a logical qubit and, by means of this 
apparatus, a breakthrough experiment shows that a universal 
set of quantum gates can be reached— Furthermore, the inter¬ 
qubit interaction can be implemented through a capacitively 
mediated dipole-dipole coupling that exploits differences be¬ 
tween charge configurations of singlet and triplet states to con¬ 
trol the interaction between logical qubits. 

In this work, we present a theoretical model to describe the 
dynamics of two singlet-triplet qubits interacting with the en¬ 
vironment. By using such a model, we are able to quantita¬ 
tively reproduce experimental results; thereby, understanding 
how singlet-triplet qubits in GaAs quantum dots interact with 
their environment. Moreover, we show that qubits are weakly 
coupled to independent dephasing channels and we reproduce 
the entanglement dissipative dynamics observed experimen¬ 
tally in Ref. i]. Such results can be used as a reference for 
further studies of open quantum system based on QDs. We 
study the entanglement dynamics as a function of the temper¬ 
ature, which plays an important role in the characterization of 
entanglement. Finally, with the correct description of the open 
quantum system, we exploit the role of the preparation time, 
i.e. the time to prepare the necessary initial superposed state, 
in the entanglement dissipative dynamics. 


The present paper is organized as follows: In Sec. II we 
present the model that describes the dynamics of the open 
quantum system based on DQDs. In Sec. Ill, we present 
the time local second-order master equation used to deter¬ 
mine the reduced dynamics of the quantum states. In Sec. IV, 
we present the concept of entanglement, measured by concur¬ 
rence together with results of our theoretical model including 
a detailed study concerning the coupling constant, the tem¬ 
perature, the environment structure, and the role of the prepa¬ 
ration time in the entanglement dynamics. Finally, Sec. V 
contains a summary of our results. 


II. THEORETICAL MODEL 

The main focus of this work is related to the study of the 
dissipative dynamics of a pair of singlet-triplet (S-Tq) qubits, 
where the information is stored in the spin state of two elec¬ 
trons. Such states can be experimentally achieved by confin¬ 
ing two electrons in each DQD system— Moreover, the two- 
level system of a logical qubit (\S) = \ t),|To) = I 4-)) can be 
isolated by applying an external magnetic field in the plane of 
the device in such a way that the Zeeman splitting makes the 
parallel spin states |r_|.) and |r_) energetically inaccessible. 

To extend such a two-level system to a two-qubit system, it 
is necessary to couple two S-Tq qubits, where the tunnelling 
between them is suppressed and their coupling is electrostatic 
(for more details, see ref. i^). Thus, the effective Hamiltonian 
for two-qubits system can be written as follows— 
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where <yx,y,z are the Pauli spin matrices, I is the identity and 
the index 1 (2) is related to first (second) qubit (hereafter, 
we use units of h = 1 ). This Hamiltonian is able to im¬ 
plement universal quantum control, which is given by two 
physically distinct local operations, x and z, and by the in¬ 
teraction between the qubits given by (g) The ex¬ 

change splitting, Ji, between [S'*) and ITq) applies rotations 
in the qubit *=1,2 around the z axis, while rotations around 
X axis are driven by a magnetic field gradient Ai?z. More¬ 
over, Ai?2 is responsible to prepare each qubit in a super¬ 
position between \S) and |To). The exchange splitting, J12, 
depends on the energy between levels of the left and the right 
DQD and it can be switched on and off during the quantum 
dynamics.- When J12 is different from zero, because of the 
Pauli exclusion principle, the \S) and ITq) states have differ¬ 
ent charge configuration, which makes the state of the first 
qubit to be conditioned to the state of the second qubit. In 
other words, when simultaneously evolving, they experience a 
dipole-dipole coupling that generates an entangled state. Fol¬ 
lowing the experimental steps,— each qubit is initialized in 
the I S') -state, then rotated by 7r/2 around the x axis when 
Ji = J12 = 0 , Ai?2 i/ 27 r Ri 30 MHz, for i=l, 2 . After this 
stage, a large exchange splitting is switched on correspond¬ 
ing to Ji/ 27 r sa 280 MHz, and J 2 / 27 r « 320 MHz. Experi¬ 
mentally, it was found that the two-qubit coupling is given by 
J12 = JiJ2^ 

To include the dissipative dynamics, we suppose that qubits 
are coupled to a bath of harmonic oscillators.- Such a cou¬ 
pling induces dephasing channels, which are the main source 
of dissipation in quantum dots-S Thus, the total Hamiltonian 
that computes the environment contribution is given by 

H = iT2-qubit + Hi, + Hint (2) 

where Hi, is the bath Hamiltonian and Hint is the interaction 
Hamiltonian. In this work, we want to determine if the envi¬ 
ronment interacts collectively or independently with the sys¬ 
tem, i.e., if both qubits noise are correlated or not. To verify 
such correlations, we suppose two distinct interaction Hamil¬ 
tonian. First, the interaction Hamiltonian is considered as a 
common dephasing environment 

Hg={aJ^ ( 3 ) 

where £ = B + with B = g^bf,. and is a complex 
coupling constant. In the common environment case, both 
qubits interact with a common bath, whose Hamiltonian is 
given by 

Hb =^^kblbk, ( 4 ) 

k 

where ujk is the frequency of the fc* normal mode of the bath 
and h\ (bk) is the creation (annihilation) operator of the reser¬ 
voir field. Superscripts of H^ and H^^ refer to collective bath. 

In the second case, the interaction Hamiltonian is composed 
of independent dephasing environments for each qubit as fol¬ 
lows 

= ( 5 ) 


where = B^^'> + B^ * with = Yl,k 9 ^k^^k 

is a complex coupling constant of the ** qubit. In this case, 

we have two independent baths described by 

«i) 

k 

where is the frequency of the fc* normal mode of the ** 

bath and (bk ) is the creation (annihilation) operator of 
the ** reservoir field. Superscripts of and refer to 
independent baths. The purpose of considering independent 
and common environment is to find a better physical model 
for the dynamics of the two DQDs open quantum system. By 
comparing our theoretical results with experimental results, 
we are able to extract such informations. 


III. MASTER EQUATION 

To calculate the dissipative dynamics, we consider a time- 
local second-order master equation, which is given by 

= - j\t'TrB{[Hj{t),[Hj{t'),pBPim} , ( 7 ) 

where pi{t) is the reduced density matrix for two DQDs, 
Hi [t) is the interaction Hamiltonian in the interaction picture, 
namely, Hi{t) = U\t)Ug{t)Hn\.UB(t)U(t), with UB{t) = 
exp (^—iHbt^ and U(t) = exp *.ff2-qubitf^ ■ Equation Q is 
valid in the regime in which the strength of the coupling, ex¬ 
pressed in frequency units, multiplied by the correlation time 
of the bath operators is much less than unity. To consider a 
collective environment in Eq.l]?]), we must employ the interac¬ 
tion Hamiltonian and the bath Hamiltonian given by Eqs. © 
and on the other hand, Eqs. © and © must be used 
for independent environments. We also suppose that the os¬ 
cillator bath density matrix ps is initially decoupled from the 
system and it is given by 


PB = — exp(-/ 3 iJ{,), 


(8) 


where Z = Tr^ 


exp(-/ 3 iTb) 


is the partition function. Here, 


P = l/{kBT), kB is Boltzmann constant, and T is the abso¬ 
lute temperature of the environment. The bath of oscillators 
is characterized by its spectral density that, in the limit where 
the number of bath normal modes per unit frequency becomes 
infinite, it can be defined asi^ 


J{uj) = r]ujexp{-uj/ujc), ( 9 ) 

where 77 is a dimensionless constant coupling that describes 
the strength of the interaction between system and environ¬ 
ment and cjJc is a cutoff frequency. 
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IV. RESULTS 

In this work, we use the quantum correlation called concur¬ 
rence to perform the analysis of our results. Concurrence is 
a well known measure of entanglement, which is broadly ac¬ 
cepted to be responsible for a set of important tasks in quan¬ 
tum information theory, such as quantum teleportationii and 
quantum key distribution!^. For two qubits, there is an analyt¬ 
ical solution to concurrencer^which is given by 

C'(p) = max{0, a/aT - \A 2 - - vAI}, (10) 

where Xi are the eigenvalues of pp listed in descending order. 
p is the time-reversed density operator, 

p= (o-W (11) 

where p* is the conjugate of p in the standard basis of two 
qubits. 

The initial state of each DQD is set to | t). then a 7 r /2 ro¬ 
tation around cc-axis is performed during the preparation time 
tpre, which put both qubits in a superposed state | t) + I i)- 
Following the experimental description given in ref. lUt], we 
use T = 50 mK, IXBzp = 2 = (tt/IO x 10^)MHz, 

Ji/2'k = 280MHz, and = 320MHz. The system 

dynamics can be obtained by numerically solving the master 
equation (Eq. (7)). We perform a systematic analysis to char¬ 
acterize the environment, i.e., we seek values for the cutoff 
frequency ujc, for the constant coupling p, and for the com¬ 
mon or the independent character of the environment that bet¬ 
ter describes the experimental data presented in ref. 0. To 
determine these characteristics, we assume a weak coupling 
between both DQDs and the environment. In such a regime, 
the dissipative process is Markovian and the cutoff frequency 
is higher than other controllable field frequencies. After de¬ 
tailed analysis (not shown here) we find, as expected, that low 
cutoff frequencies are unable to reproduce the experimental 
data and we checked that all results presented in this work do 
not change significantly if Wc > 2 x lO'^MHz. Thus, the cutoff 
frequency Wc = 2 x lO^MHz is fixed hereafter. 

We begin our analysis by comparing the dynamics for 
the difference of the descending sorted eigenvalues (DDSE) 
Ai — A 2 — A 3 — A 4 of the matrix pp, which is equal to the con¬ 
currence C{p) when it assumes positive values. In Eig. ([T]i, 
we plot Ai — A 2 — A 3 — A 4 as function of time, considering 
a common environment coupled to both qubits and assuming 
different values for the coupling constant p, for the experi¬ 
mental temperature T = 50 mK. One can see in Eig. ([T]i that 
the DDSE oscillates on time and it has an envelope function 
which decays faster as the constant coupling is increased. We 
also perform an analysis of DDSE as a function of time for 
a fixed coupling constant p = 3xl0“® and different tem¬ 
peratures, which is shown in Eig. The increasing of the 
temperature has a similar effect as the increasing of the con¬ 
stant coupling (Eig. O; i-e., the higher the temperature, the 
faster the decay of the envelope function as a function of time. 
Based on results of Eig. ([TJ and Eig. (|2]), we conclude that the 
description of common environment for both DQDs does not 
describe the experimental results observed in ref. JH] because 



FIG. 1: (Color online) Numerical solutions for the DDSE of the ma¬ 
trix pp as a function of time, for a fixed temperature T = 50 mK and 
a fixed cutoff ujc = 2 x lO'^MHz, considering a common environ¬ 
ment coupled to both qubits. The dot-dashed (blue) line represents 
the DDSE for 77 = 10“®, while the dotted (orange), the solid (black), 
the double dotted (red), and the double dot-dashed (magenta) respec¬ 
tively correspond to 77 = 2 x 10~®, 77 = 3 x 10~®, 77 = 4 x 10~®, 
and 77 = 5 x 10“®. 



FIG. 2: (Color online) Numerical solutions for the DDSE of the ma¬ 
trix pp as a function of time, for a constant coupling 77 = 3 x 10“® 
and a fixed cutoff ujc = 2 x 10^MHz, considering a common en¬ 
vironment coupled to both qubits. The dot-dashed (blue) line repre¬ 
sents the DDSE for T = 0 mK, while the dotted (orange), the solid 
(black), the double dotted (red), and the double dot-dashed (magenta) 
respectively correspond to T = 10 mK, T = 50 mK, T = 100 mK, 
and T = 200 mK. 


the DDSE do not present negative values, as found experimen¬ 
tally. Eurthermore, the concurrence (positive values of DDSE) 
theoretically obtained presents a sudden-birth!^ which is not 
observed in the experiment. 

To provide a better description of such a experimental re¬ 
sult, we also assume independent environments for each qubit. 
We choose the experimental temperature T = 50 mK, in 
Eig. m, and we plot the DDSE as function of time, con¬ 
sidering different coupling constants. Remarkably, for 77 = 
3 X 10“®, tUc > 2 X 10"*^MHz, there is a very good match¬ 
ing between experimental results (open circles in Eig. Q) 
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FIG. 3: (Color online) Numerical solutions for the DDSE of the ma¬ 
trix pp considering independent environments for each qubit, as a 
function of time for a fixed temperature T = 50 mK and a fixed cut¬ 
off LUc = 2 X 10'*MHz. The dot-dashed (blue) line represents the 
DDSE for p = 10“®, while the dotted (orange), the solid (black), 
the double dotted (red), and the double dot-dashed (magenta) re¬ 
spectively represent the results for 77 = 2 x 10~®, p — 3 x 10~®, 
77 = 4 X 10“®, and 77 = 5 x 10“®. Open circles denote the DDSE 
extracted from experimental datai^ 



EIG. 4: (Color online) Numerical solutions for the DDSE of the ma¬ 
trix pp as a function of time, for a constant coupling p — 3 x 10~® 
and a fixed cutoff uic = 2 x lO^MHz, considering independent en¬ 
vironments for each qubit. The dot-dashed (blue) line represents the 
DDSE for T = 0 mK, while the dotted (orange), the solid (black), 
the double dotted (red), and the double dot-dashed (magenta) respec¬ 
tively represent the results for T — 10 mK, T — 50 mK, T = 100 
mK, and T = 200 mK. 


and DDSE extracted from the reduced dynamics obtained for 
independent environments for each qubit. The experimental 
data for the DDSE, plotted in Eig. Q, achieves its largest 
value DDSEr:; 0.45 for t = 150 ns; also, concurrence (pos¬ 
itive values of the DDSE) is null for t > 250 ns. The time 
where the concurrence is maximum can be related to term 
® cri^^ that couples both qubits, thus its amplitude 
J 12 is connected to the time where the maximum entangle¬ 
ment (concurrence) occurs. Such relation is J 12 = 7r/re„7, 


EIG. 5: Numerical solutions for the maximum value of DDSE, con¬ 
sidering independent environments, as a function of the amplitude of 
R — Ji 2 /Ti 2 ^ and for different preparation times tprep- The solid 
(black) is the result for tprep = 4 ms, while the dotted (orange), the 
double-dotted (red), and the double dotted-traced respectively show 
the results for tprep = 8 ms, tprep = 32 ms, and tprep = 64 ms. In 
the inset, we plot the DDSE of the matrix pp as a function of time 
following the same preparation times and curve description used in 
the main panel. 



FIG. 6 : Numerical solutions for the maximum value of DDSE as a 
function of tprep for T = 50 mK, 77 = 3 x 10“®, cOe = 2x 10"* MHz, 
R = 50, and independent environments assumption. 


which gives Tent = 150n8^ in such a case. 

By means of the correct description of the environment, we 
can analyze the entanglement dissipative dynamics through 
our theoretical model. We begin our analysis focusing on the 
role of the temperature in the dissipative process. In Fig. 
we fixed the parameters as the ones that better fit the exper¬ 
imental results and plot the DDSE as a function of time, for 
different temperatures. For T = 0 K, we see that the electro¬ 
magnetic vacuum has a negligible effect on the deterioration 
of the entanglement, which only oscillates as a function of 
time. Moreover, from Fig. we see that a change of just 10 
mK significantly influences the entanglement dynamics and 
for T = 200 mK, the entanglement definitely disappears. 

We also study the role of the preparation time tprep, i-e. the 
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time to prepare the superposed state. In the experiment, each 
qubit is initially prepared in the | t)-state, then rotated by 7r/2 
around the x-axis when Ji = J 12 = 0 to achieve the super¬ 
posed state I t) + I i)- This rotation around x-axis is driven 
by a magnetic field gradient During the preparation 

time, the system is interacting with the environment, thereby 
affecting the entanglement dynamics. To probe such an effect, 
we study the role of tprep in the maximum value of entangle¬ 
ment during the dissipative dynamics. In the experimental re¬ 
sult, for example, we have tprep ~ 8 ns, ~ 7r/150ns, 
and the maximum obtained entanglement is around 0.45 for 
T = 50 mK. Naturally, another crucial aspect to maximize 
the entanglement between qubits is the exchange coupling be¬ 
tween each DQD. The exchange coupling J 12 , which can be 
increased by controlling the dipole-dipole interaction, deter¬ 
mines the time for achieving the maximally entangled state 
and, the faster the maximally entangled state is prepared, the 
smaller the environment perturbation. 

To illustrate the role of the preparation time in the entangle¬ 
ment dynamics, we plot the maximum value of entanglement 
between the qubits as a function of R = J\ 2 l for a dif¬ 
ferent preparation time tprep in Fig- ©. As expected, such 
results show an enhancement of the entanglement when J 12 
is increased and when tprep is decreased. This behavior is re¬ 
lated to the fact that a maximally entangled configuration is 
faster achieved for larger Jyi but the interaction with the envi¬ 
ronment disturbs the ideal initialization during the preparation 
time tprep- To understand how the preparation time affects the 
maximally entangled state, in Fig. we plot the maximum 
value of concurrence as a function of tprep for a fixed R = 50. 
One can notice that the maximum concurrence linearly de¬ 


creases as a function of tprep- Such results show what are the 
parameters to enhance the entanglement in S'-Tq qubits. 


V. CONCLUSION 

In summary, we proposed a model to describe a open quan¬ 
tum system composed of two DQDs. By employing such a 
model, we were able to compare the case where both DQDs 
are coupled to a common environment to a situation where 
each DQD is coupled its own bath of harmonic oscillators. By 
performing a systematic analysis of the environment descrip¬ 
tion, we found that the independent environment case agrees 
with the experimental data shown in ref. ||3l- Moreover, we 
study the role of the temperature and of the preparation time 
in the dissipative dynamics. We show that the temperature 
presents a crucial role in the entanglement evolution. Fur¬ 
thermore, we show that the maximally entangled state decays 
linearly with the preparation time. 
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